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Abstract. We derive general properties of the finite-size scaling of probability density functions and show 
that when the apparent exponent f of a probability density is less than 1, the associated finite-size scaling 
ansatz has a scaling exponent r equal to 1, provided that the fraction of events in the universal scaling 
part of the probability density function is non- vanishing in the thermodynamic limit. We find the general 
result that r > 1 and r > f . Moreover, we show that if the scaling function G(x) approaches a non- zero 
constant for small arguments, Yim x ^oG{x) > 0, then t — f. However, if the scaling function vanishes 
for small arguments, lim^-^o G{x) = 0, then r = 1, again assuming a non-vanishing fraction of universal 
events. Finally, we apply the formalism developed to examples from the literature, including some where 
misunderstandings of the theory of scaling have led to erroneous conclusions. 
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Power-law probability densities are pervasive in the lit- 
erature on critical and scale- invariant systems [1-3]. The 
supposedly sole scaling parameter in a finite system at the 
critical point is the upper cutoff s c , diverging in the ther- 
modynamic limit, L — > oo, and the probability density 
function (PDF) P(s; s c ) for an observable (an event size 
or the order parameter) s £ [0, oo[ obeys simple finite-size 
scaling (FSS), that is, 



P(s;s c ) = as T Q{s/s c ) for s,s c >s , 



(1) 



where so is a constant lower cutoff of the PDF. Equa- 
tion (1) is valid only in the region s 3> sq where the 
corrections to scaling are negligible [4]. The dimension- 
ful parameter a is a so-called non-universal metric factor 
[5], t is a universal (critical) scaling exponent, Q is the 
universal scaling function that decays sufficiently fast for 
s ^> s c to ensure that, as one would expect, all moments 
of the PDF are finite for finite system size, L < oo. The 
upper cutoff s c is the characteristic size of s in a finite 
system. Usually s c = bL D to leading order, where b is an- 
other non-universal metric factor and D is the universal 
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spatial dimensionality of the observable s. Again, in gen- 
eral, there are sub-leading orders, that is, corrections to 

scaling of the form s c = bL D (1 + 6iL" Wl + b 2 L~ u ' 2 H ) 

with < u>i < LU2 < • • •, which can be safely ignored in 
the analysis of the asymptotes. The scaling exponent r 
is uniquely defined by Eq. (1) and is associated with the 
power-law decay of the distinctive onset of the rapid de- 
cay P(s c ; s c ) — aQ(l)s c ~ T as a function of s c . If data are 
consistent with the FSS ansatz (1), we can perform a data 
collapse: By plotting the transformed PDF s T P(s; s c ) vs. 
the rescaled observable s/s c , all data for s,s c 3> sq col- 
lapse onto the same curve representing the scaling function 
Q, see Figure 1. However, due to the presence of the metric 
factor a, the scaling function Q is determined by Eq. (1) 
only up to a prefactor. This ambiguity can be resolved, 
for example by fixing 0(1) = 1 [6]. 

The upper and lower cutoffs define a scaling region 
sq <C s <C s c where the PDF shows essentially power-law 
behaviour, rather informally 



P(s; s c ) oc s T for sq <C s <C s c 



(2) 



a Present address: Department of Mathematics, Imperial 
College London, 180 Queen's Gate, London SW7 2AZ, United 
Kingdom 



which means that s~ r is the leading order of P{s; s c ) in the 
region where s > so coexists with s <C s c . The apparent 
exponent f is the slope of a straight-line fit to the PDF 
data in the scaling region sq -c s <C s c when plotting 
logP(s;s c ) vs. logs, see Figure 1(a). 
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Let us reflect upon the FSS ansatz (1) and the infor- 
mal Eq. (2). Equation (2) is informal for various reasons. 
First, the "oc" suggests a proportionality without specify- 
ing what parameters a possible prefactor would be allowed 
to depend on. Second, the upper cutoff appears in the con- 
dition s<s c , but s c does not enter Eq. (2) explicitly as it 
does in Eq. (1) through its role as a characteristic scale in 
the dimensionlcss argument of the scaling function Q. Fi- 
nally, and most importantly, the scaling function Q is not 
present. Hence, to derive Eq. (2) from Eq. (1), one has to 
make assumptions about the behaviour of Q in the scal- 
ing region sq <C s <C s c . If the scaling function is "almost 
constant", then r = f . However, if the scaling function is 
a power law itself, then r ^ f . To illustrate that, let us 
define a cutoff function G, that, by definition, converges 
to a non-vanishing value for small arguments and decays 
rapidly for large arguments, and assume that the scaling 
function is of the form Q{x) = x a Q(x). Since the cutoff 
function Q is constant and non-zero in the scaling region, 
the PDF behaves like as~ T (s/s c ) a Q(0) for s <C s < s c , 
so that the scaling exponent r in Eq. (1) is related to the 
apparent exponent f in Eq. (2) via r — a = f. 

One of the key problems found in the literature is the 
identification of the scaling exponent r with the apparent 
exponent f; however, they differ if the scaling function has 
a power-law dependence on its argument. First, we illus- 
trate the consequences of such a scaling function further, 
concluding that if the apparent exponent f in Eq. (2) is 
less than unity, then the scaling exponent r in Eq. (1) is 
exactly 1. Next, this result is derived in general by putting 
the formalism, Eqs. (1) and (2), on a sound mathematical 
basis, generating further general properties in the process 
such as t > 1, t > f and f > 1 r = f. Finally, 
we review examples from the literature that illustrate the 
formalism developed, including some that suffer from a 
confusion of r and f. 



1.1 Example 
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Fig. 1. (a) The PDF P(s;s c ) vs. s with an apparent expo- 
nent t — 0.8 (indicated by a straight line with slope —0.8). 
Since f < 1, the prefactor a decreases with increasing cutoff, 
shifting the graphs downward the ordinate, see Eq. (4b). (b) 
Data collapse. The transformed PDF sP(s; s c ) vs. the rescaled 
observable s/s c - For s, s c S> so the curves collapse onto the 
graph for the scaling function G(x) — x 1 ~ T Q(x) where the cut- 
off function Q{x) — (1 — x) r 9(l — x), see Eq. (3). For small 
arguments the scaling function behaves like a power law x ' 2 , 
as indicated by a straight line with slope +0.2. Example taken 
from Ref. [7]. 



To demonstrate the first result explicitly, we assume that 
the PDF for the observable s e [s , oo[ has the form 

P(s;s c ) = as~ T <3 (s/s c ) for s > s = 1, (3) 

see Fig. 1(a). Motivated by a previous study [7], in this 
example, the cutoff function is assumed to be Q{x) = 
(1 — x) T 6(l — x) with 8 being the Heaviside step function. 
(A simpler, but more artificial example would be to con- 
sider a cutoff function identical to a pure Heaviside step 
function.) Along the lines of the considerations above, this 
cutoff function converges to unity for small arguments and 
vanishes for all x > 1. The apparent exponent, the slope of 
P(s; s c ) in the scaling region as shown in Fig. 1(a), there- 
fore is the exponent f in Eq. (3). A naive comparison of 
Eq. (I) and Eq. (3) suggests that the scaling exponent r is 
the slope of a straight-line fit to P(s; s c ) vs. s in a double 
logarithmic plot, but this is incorrect if f < 1, as we will 
see in the following. 



In this particular example (3), in the limit of large s c 
the scaling form of the PDF is expected to account for all 
s, so that in this limit 

/>OC 

1= / dsP(s;s c ) 

Js = l 

= (a-± I forf>l, (a) 

\as c 1 - f r(i-f)r(i + f)forf <i. (b) y i 

For simplicity, we ignore the case f — 1 which contains 
logarithmic corrections. Equation (4) determines a, which 
for f < I cannot be a constant but has to vanish like 
s c r_1 in the limit of large s c . Substituting a into Eq. (3) 
and rearranging, the PDF reads 

P(s;s c ) (5) 

= f (f-l)s- f g(s/s c ) forf>l, 

~ I ra-f/ni+f) 5 ' 1 ( s / s c) 1_T Q (s/s c ) for f < 1. 
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By comparing Eq. (5) with the FSS ansatz (1), we see 
that when the apparent exponent f > 1, the metric fac- 
tor a = (f — 1), the scaling exponent r = f, and the 
scaling function Q(x) = Q(x). However, if the measured 
apparent exponent f < 1, the metric factor a = (T(l — 
f)r(l + f)) -1 , the scaling exponent r = 1, and the scal- 
ing function Q(x) = x 1 ~ T G(x). For the FSS ansatz (1), 
it is a necessity that the metric factor is asymptotically 
independent of s c . Dropping this constraint would render 
the definition of the scaling exponent r meaningless: One 
may rewrite as~ T Q(s/ s c ) = as" s~( T+a \s / s c ) a Q(s / s c ) for 
any a. The factor (s/s c ) a could be absorbed into the scal- 
ing function. Were it allowed to absorb the power s" in 
the metric factor, it would be impossible to differentiate 
between r and r + a. 



2 Derivation 

Now we present a more general and rigorous derivation of 
the result illustrated above. 



2.1 Definitions and notation 

The FSS ansatz (1) can be recast more formally as follows: 
For all e > there exists a constant lower cutoff so and 
an S so that for fixed n > 0, the relative error of the 
nth moment from the approximate PDF of the observable 

s e [0, oo[ 

p , s , = I as~ T g(s/s c ) for s > s , (6a) 
\ f(s, s c ) otherwise (6b) 

compared to the exact nth moment of the PDF approxi- 
mated by -P app (s; s c ) is less than or equal to e for all s c > S. 
Here, we have introduced /(s,s c ), the non-universal part 
of the PDF. Both functions, /(s,s c ) and G{x), are non- 
negative, and the latter is assumed to be non-zero at least 
in a finite range of x. We assume s > 0, but the case 
so = is briefly discussed as well. 
The exact nth moment is given by 

(s n ) = / n (*o, *c) + a S l +n - T g n ( So / Sc ) + 0(e) (7) 

where 

r s o 

fn(so,s c )= / dss n f(s,s c ), (8a) 
Jo 

P CO 

9n(s /s c ) = dx x n - T G(x). (8b) 

J sq/s c 

Both /„ and g n are non-negative and g n is a monotonically 
non- increasing function of so/s c . It is worth noting that 
the range over which the scaling region is probed widens 
as s /s c approaches 0. 



2.2 Scaling exponent r > 1 

Since go does not decrease with s c , the normalisation con- 
dition 

1 = (s°) = /o(*o, 5 C ) + as c 1 - T .g (s /s c ) + 0(e) (9) 

implies r > 1, otherwise asl~ T go(so/s c ) would diverge as 
s c — > oo. Therefore, the scaling exponent r > 1. 



2.3 If Q is continuous, the limit Q = km x ^oG{x) 
exists and is finite 

The behaviour of the scaling function Q in the limit of 
small arguments s/s c enters into the scaling of the mo- 
ments in the asymptotic limit of large upper cutoff (the 
thermodynamic limit). We will assume that G(x) is con- 
tinuous within a finite range ]0, q] with q > and that 
the integrals in Eq. (8) exist. Without these (reasonable) 
assumptions, the discussion of the behaviour of the inte- 
gral g n becomes too complicated. The continuity implies 
further that G(x) cannot diverge with small arguments: 
If it did, go(so/s c ) diverges faster than s£ _1 for r > 1 
and at least logarithmically fast for r = 1, contradicting 
the normalisation condition Eq. (9). Therefore, if G{x) 
is continuous, it cannot be divergent in small arguments. 
Together with the continuity, the limit Go = lim^o G{x) 
therefore exists and is finite. 



2.4 If t = 1 then Go = 

We consider the case Go > 0: If G{x) is non-zero and 
bounded from above for small x > 0, the integral g n (so / s c ) 
is logarithmically divergent in s c when I + n — r = 0, 
has leading order s c ~' 1+ " _T ' ) when 1 + n — r < and 
converges to a non-zero value with increasing s c when 
I + n — t > 0. The normalisation condition Eq. (9) im- 
plies asl~ T go(so/s c ) < 1, so that Go > violates the nor- 
malisation condition if r < 1. For r > 1, the divergence 
s c ^ ^ of go is cancelled by the prefactor s]r T . Therefore, 
if Go > then r > 1. Negating this result yields: If t = 1 
then <?o = 0. 



2.5 If the lower cutoff s = then r = 1 and Go = 

The result above can be amended by the observation that 
if the lower cutoff s = then all non-universal /„ vanish 
and the normalisation condition Eq. (9) reduces to I = 
asl~ T go(0) + 0(e), which can hold only if r = 1. Moreover, 
So(0) — J °cla; x~ 1 G{x) can converge for continuous G(x) 
only if Go — 0, consistent with our findings in Section 2.4 
Therefore, if the lower cutoff vanishes, so — 0, then r = 1 
and Go = 0. Generally, the converse does not hold, that is, 
Go = does not imply s = 0. However, Go > or r > I 
implies s > 0. 
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2.6 If So = then r = 1 

We consider the case Go — 0: Here, we focus on g . If r > 1, 
then g (so/s c ) diverges slower than (so/s c ) 1 ^ T in decreas- 
ing so/s c . This can be proven by assuming the opposite, 
implying that there exists a c > so that J s / s dx x~ T Q{x) > 

c (sq/sc) 1 t for all sq/s c below a threshold determined by 
c. The right hand side of this relation can be written as 
(r — 1) J s / s dx x~ T . Independent of its prefactor, for suf- 
ficiently small x the integrand represents an upper bound, 
given that Q(x) is less than any 5 for sufficiently small x. It 
is easily shown that the above relation therefore eventually 
breaks down, implying that the product asl~ T g (so/s c ), 
the fraction of events within the universal scaling part of 
the PDF P(s;s c ), vanishes in increasing s c if r > 1. No 
such argument exists if r = 1. Mathematically, a van- 
ishing fraction of such "universal events" constraints nei- 
ther G(x) nor r. However, it makes sense to demand some 
weight in the scaling part of the PDF on physical grounds. 
Therefore, if Go vanishes and the fraction of events in the 
scaling part of the PDF is finite in the thermodynamic 
limit, then the scaling exponent r = 1 . If Go vanishes and 
the scaling exponent r > 1, then the fraction of events in 
the scaling part of the PDF vanishes in the thermodynamic 
limit. Hence, we are faced with the almost paradoxical sit- 
uation that all moments with n > r — 1 are determined 
by the scaling part of the PDF, the measure of which, 
however, vanishes in thermodynamic limit. In this case, 
the normalisation condition Eq. (9) entails that /o(so7 s c) 
converges to 1 in the limit of large s c , that is, the lower 
cutoff so cannot vanish. Finally, we note that, if Go = 0, 
then g n (so/s c ) cannot diverge faster for sq/s c — ► than 
in the case Go > since G(x) does not diverge for small 
arguments. Therefore, if Go = then g n (s o / s c ) converges 
for all n > r — 1 . 



2.7 If Go > then r = f 



The apparent exponent f is formally defined as 

lim lim P ^ S; Sc } \ f = 1 for all A > 0, (10) 

s^oc sc^oc P(S; S c ) 

if the limits and such an exponent exist, which is, for ex- 
ample, not the case if G (x) converges to faster than any 
power law. The limit s — > oo guarantees that the limit 
s c — > oo is taken in the scaling region, so that P(s; s c ) can 
be replaced by the scaling part Eq. (6a): 



1 for all A > 0. 



(11) 



If Go > 0, the fraction in Eq. (11) converges to 1 and 
for the left hand side to be equal to unity for all A, the 
exponent f — r vanishes. Therefore, if Go > then r = f 
or, equivalently, t ^ f => Go = 0. 



2.8 Generally r > f , if f < 1 then r = 1, and if f > 1 
then r = f 

Equation (11) implies that £(As/s c ) scales like A T_T and 
therefore diverges in small arguments if r < f. However, 
as noted previously, G cannot diverge for small arguments. 
Therefore, we have the general result r > f. 

Moreover, the cutoff function G(x) = x T ~ T G(x) can 
then be expected to converge to a (non-vanishing) con- 
stant as x — > 0, and imposing that, one can define uniquely 
the apparent exponent f through 

P(s;s c ) = a( Sc ) s - f G( s /s c ) fors>s - (12) 

Using the above two arguments, the scaling exponent 
t and the apparent exponent f can be related as follows: 
If f < 1 then t/t since r > 1, implying Go = 0. As 
shown above, Go = implies r = 1 assuming that the 
scaling part of the PDF is non-empty. Therefore, if f < 1 
then t = 1. 

Now, assume to the contrary that t ^ f > 1. Hence 
Go = and again, assuming a non-empty scaling part of 
the PDF, t = 1. If f = 1 that clashes with r ^ f, if f > 1 
that clashes with r > f . Therefore, iff > 1 then t = f, 
unless the total measure of the scaling part of the PDF 
vanishes asymptotically. 



3 Discussion 

In this section, we first discuss some implications of the 
above derivations and then present some examples from 
the literature, where a more thorough scaling analysis 
leads to revised results. 

In many models of equilibrium [8] as well as non-equi- 
librium [9] statistical mechanics, within the FSS regime, 
the second moment (s 2 ) of the total order parameter 1 
scales like its first moment squared (s) 2 , which is usually 
linked to so-called hyper-scaling. This implies that their 
ratio converges in the thermodynamic limit, that is, 



lim 



(f!> 



(13) 



where r is a finite number. We now show that such fi- 
nite r implies r = 1, assuming that Eq. (1) holds in these 
models, by showing that if r > 1, then the ratio diverges 
asymptotically. To invoke the argument "if t > 1 then 
Go > 0", we need to impose that a non- vanishing frac- 
tion of the total order parameter distribution is governed 
by finite-size scaling. In equilibrium statistical mechanics, 
that means that the singular part of the partition sum 
does not vanish asymptotically. Moreover, we assume that 
both the first and the second moment diverge in the ther- 
modynamic limit. This is important, because r may be 



1 We distinguish between the total order parameter s and the 
order parameter density s/L d . The analysis is identical for the 
order parameter density, however, the upper cutoff then van- 
ishes asymptotically which would require a revision of previous 
arguments. 
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finite even for Go > simply by convergence of the first 
two moments if r > 3, as can be shown using Eq. (7) and 
the divergence of g n as discussed in Sec. 2. 7. Based on the 
arguments presented there, if the first moment diverges, 
then r < 2. Using the same arguments, one can now show 
that for 2 > r > 1 the ratio r diverges like s^ _1 and for 
r = 2 it diverges like s,T _1 with logarithmic corrections. 
Thus asymptotically finite r is incompatible with t > 1, 
if the first moment diverges and the fraction of the total 
order parameter distribution that is governed by the scal- 
ing part is finite. The only alternative is t = 1, in which 
case Go = and therefore all g n are finite for all n and 
r = u~ l g2/9i horn Eq. (7). The case r = 1 therefore is 
anything but special, describing most of standard criti- 
cal phenomena. It is worth noting that r = 1 implies a 
dimensionlcss metric factor a. 

We illustrate the observation above by considering the 
finite-size scaling of the cluster number density 2 in per- 
colation [2] (i.e., the number of clusters of size s per site 
in a system of linear size L) 

n{s; S r°) = a perc s- rPer > erc (s/s^) , (14) 

where sg orc = 6 porc L D , D = \j(ya) with v and a being the 
critical exponents describing the divergence of the correla- 
tion length and the characteristic cluster size, respectively, 
and 6 porc being a metric factor. The exponent r porc is differ- 
ent from 1, known to be exactly 187/91 in two-dimensional 
percolation. However, n(s;sg orc ) is not the distribution of 
the total order parameter, which is the number of occu- 
pied sites in the largest cluster. Multiplying n(s:s porc ) by 
L d , one obtains the average number of clusters of size s 
in a system with L d sites. The probability of the largest 
cluster being of a particular size must be contained in 
this distribution: If n(s; s porc ) was the site- normalised den- 
sity of the largest cluster only, then L d n(s; s porc ) was the 
probability that the largest cluster has size s. If n(s; s porc ) 
is the site-normalised histogram of all clusters exclud- 
ing the largest, one could therefore write L d n(s; s porc ) = 
L d n(s; s porc ) + P(s; s porc ) with P(s; s porc ) being the distri- 
bution of the total order parameter. One might speculate 
that P(s;s porc ) is responsible for the characteristic bump 
in the histogram n(s; s porc ). Assuming that P porc (s; s perc ) 
follows the same scaling as L d n(s; s porc ) itself, one has 

ppe r c (s; = Ld s -r—g P ^ {s/sr ^ (15) 

with s porc = b poTC L D . Absorbing L d into a redefinition of 
the scaling function G v '" c {x) = x^ d / D G po "{x) the distri- 
bution of the order parameter scales like 

ppc rc(s . s pe re) = a ^- T P- gperc(s/sro)) (lg) 

where the exponent of s turns out to be d/D — t p<! '° = — 1 
by standard scaling laws [2], indeed as expected from the 
previous analysis. 

2 In the following, we consider finite systems at occupation 
probability p tuned to the critical value p c , so that we do not 
have to exclude from the cluster number density the largest, 
infinite or percolating cluster(s). 



3.1 Examples 

Finally, we want to present some examples from the liter- 
ature, where a misunderstanding of the theory of scaling 
has lead to some confusion or even erroneous conclusions. 
In [10] the cluster number density in steady state two- 
phase flow in porous media is investigated. As discussed 
above, in percolation the cluster number density n(s) is 
the site normalised density of clusters of size s, so that 
sn(s) is the probability that a randomly chosen site be- 
longs to a cluster of size s. The histogram is expected to 
follow scaling, n(s) = a pcrc s rP ° rc 5 pcrc (s/s pcrc ), see Eq. (14). 
The sum p = J2 S sn i s ) is the fraction of sites belonging to 
any cluster, i.e. the fraction of occupied sites. This sum 
is bounded and non-zero in non-trivial cases even if the 
largest clusters are excluded from n(s), so that the re- 
sults for P(s; s c ) derived in this paper, translate into cor- 
responding results for sn(s). Performing the sum suggests 
p cx (s pcrc ) 2 ~ rPcrc , which means that p vanishes asymptot- 
ically for t p °" > 2. This, however, is a misunderstanding 
of the theory of scaling, because it ignores the existence 
of a finite lower cutoff Sg orc . In fact, r pcrc > 2, correspond- 
ing to t > 1 as derived earlier, is a necessary condition 
for p to be finite. The numerical finding of r porc w 1.92 in 
[10] for a cluster number density therefore is theoretically 
impossible - in fact it is based on a straight line fit in a 
double logarithmic plot and therefore it is f and not r 
that is estimated. 

The term self-organised criticality refers to the ten- 
dency of slowly driven non-equilibrium systems with many 
degrees of freedom to spontaneously reach a critical state 
where the PDF of relaxational event sizes obey FSS. There 
are only a few controlled laboratory experiments on scale- 
invariance in slowly-driven non-equilibrium systems. In 
one such experiment, the avalanche-size PDF was mea- 
sured in slowly driven piles of rice [11]. For piles with elon- 
gated grains, the avalanche-size PDF obeys FSS. However, 
it was claimed that a system with more spherical grains 
is non-critical, since the PDF of the avalanche size s in a 
system of size L is consistent with P(s;L) cx L~ 1 G(s/ L) 
where the function G(x) oc cxp [— (x/x*) 1 ] with 7 w 0.43 
and x* w 0.45 was identified as a scaling function rather 
than a cutoff function. This is a misunderstanding of the 
theory of scaling. The PDF can be recast into the FSS 
form P(s;L) cx s~ 1 G(s/L) where the scaling function 
G(x) = xG(x) increases linearly with x for small argu- 
ments, decaying (stretched) exponentially fast for large 
arguments. Hence, also for the piles with more spherical 
grains does the avalanche-size PDF obey FSS but with a 
scaling exponent r = 1 and Go = 0. Therefore, the occur- 
rence of self-organised criticality does not depend on the 
details of the system as wrongly claimed in [11]. 

Complex networks have attracted intense interest from 
the statistical physics community in particular because 
they often exhibit scale invariance. It has been suggested 
that so-called duplication-divergence- mutation models, where 
nodes correspond to proteins and links correspond to pair- 
wise interactions between proteins, are good candidates 
to describe the evolution and large-scale topological fea- 
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tures of real protein-protein interaction networks. A recent 
study [12] derives asymptotic properties for such models 
and it is reported that the degree PDF for a node to have 
degree k at time t is P(k;t) oc t- 1 (k/t) k ^- 1 G(k/t) where 
fc min is the lower non-zero degree in the initial configura- 
tion, Q a cutoff function, and it is claimed that the scaling 
exponent is fc min — 1 [12]. However, this is a misunder- 
standing of the theory of scaling. A careful analysis reveals 
that the PDF can be recast into the FSS form P(k; t) = 
ak^ 1 Q(k/t) where the metric factor a — fc min (l + fc min ) and 
the scaling function Q(x) — x fcmin (l — x)9(l — x) increases 
like a power law with exponent fc min with x for small argu- 
ments, decaying fast for large arguments. Hence, the FSS 
scaling exponent r = 1, independent of fc min and Go = 0. 

Earthquake statistics has been a prominent subject of 
statistical analysis by means of scaling arguments (e.g. 
[13]). In [14], the distribution P m ^(Ar) of distances of 
successive earthquakes with a magnitude greater than m 
within a cell of linear L is considered. The analysis shows 
that P m . L {Ar) = L~ 1 f(Ar/L), where f(x) oc a;" - 6 . As 
that implies P m ^(Ar) oc Ar~ oe , the distribution seems 
to be asymptotically non-normalisable. This conclusion, 
however, is a misunderstanding of the theory of scaling 
and a confusion of f and r: The PDF P m ^(Ar) was non- 
normalisable only if t < 1, but in fact r = 1, which can be 
seen by rewriting the scaling of the PDF as P m ^(Ar) = 
Ar~ 1 Q(Ar/L) with Q(x) = xf(x). The scaling function 
Q(x) converges to for small arguments like a; 1 " 6 , as 
expected from the discussion above. In the introduction 
we showed that Q(x) oc x a implies f = r — a, so that in 
fact t — 0.6 in the present case. 

In summary, given the probability density function 
P(s; s c ) obeys the FSS ansatz (1), if the apparent expo- 
nent f of P(s; s c ) in Eq. (2) is less than unity, then the 
scaling exponent r is unity and the scaling function van- 
ishes for small arguments, lim^^o G( x ) = 0. If f > 1, then 
the scaling exponent equals the apparent exponent, r = f . 
If r > 1, then the scaling function approaches a non-zero 
value for small arguments, lim^o Q(x) > 0. 

GP thanks Mervyn Freeman, Nicholas Watkins and Max Werner 
for useful discussions. 
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